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I. INTRODUCTION 


The molecular species C0Qz and H20O are major atmospheric 
constituents which critically control the infrared radiative 
transfer in the earth's atmosphere. Their absorption bands 
are distributed over a wide spectral range, extending from 
the near to the far infrared region of the electromagnetic 
spectrum. These species have been extensively studied since 
infrared spectroscopy technology came into existence; however, 
there are still gaps in our knowledge of the infrared absorption 
of these two molecules. In particular, our understanding of 
their absorption at high temperatures remains incomplete. 

The bands observable at room temperature are restricted to 
transitions which originate from a lower state of a relatively 
small excitation energy—typically from a vibrational-rotational 
energy below 2000 cm-l. Transitions at room temperature are 
weak since these lower states are insufficiently populated. 

They cannot be studied well even with measurements through 


extremely long absorption paths. These are the transitions 


which originate from lower states of high vibrational-rotational 


energy and become observable when the gas temperature is raised 

to increase their population density. These transitions, 4 

generally referred to as the hot bands, are poorly known : 

because of the paucity of measurements. i" 
Various experimental difficulties are encountered in the | 


spectroscopic measurement of gases at elevated temperatures. 


The original AFGL compilation of the atmospheric line parameters 


was assembled primarily for application to room temperature 


calculations. The basic line parameters were derived from 


Consequently, the accuracy which the 


room temperature data. 


data provided was generally unsuitable for spectral synthesis 


calculation for high temperatures. The line parameters 


compiled for the hot bands of these molecules do not provide 


an accuracy sufficient for line position extrapolation to 


high temperature. The present experimental effort was begun 


with the intent to solve this difficulty by measuring the 


absorption of these molecules at elevated temperatures and 


obtaining improved line parameters which could yield more 


accurate results in spectral syntheses at high temperatures 


and high resolution. 


The excitation which occurs in a heated gas maintains 


thermal equilibrium of the molecular system; the distribution, 


or population as a function of energy level, is controlled 


by the Boltzmann factor e~E/KT, where E is the vibrational- 


rotational energy of the state in question. There would be 


no preferential excitations among various vibrational levels. 


If two states which belong to two different vibrational states 


have a similar energy, E) ~ Ej, both states are excited at 


an equal rate. A marked increase in the observable transitions 


consequently occurs along two directions as the temperature 


of the molecular system is increased; there is a marked 


increase in observable vibrational transitions as well as 


in rotational transitions within a vibrational transition. 


The spectral structure consequently increases in complexity 
as the gas temperature is raised. 

The increased complexity in the spectra poses two 
problems: one is a requirement for increased spectral resolu- 
tion, and another is for the identification of the individual 
spectral transitions. The spectra must be observed with 
adequate spectral resolution for the isolation of most of 
the transitions of interest. The spectral coverage must be 
broadened since more rotational lines are excited within a 
vibrational band. The problem of spectrometry is, for our 
case, adequately solved by use of the technique of Fourier 
spectroscopy. The second problem concerning identification 
of the observed spectral transitions was found more difficult 
than expected. To overcome this difficulty the spectral 
data were taken at temperature steps of 200K, i.e., at 600K 
and at 800K. We hoped to follow the excitation of the hot 
bands in the data taken at these temperatures. For the H)0 
data, the temperature step taken was adequate, while it was 
found inadequate for the CQ) data. The increase in observable 
CO2 lines from 600K to 800K was overwhelming; at 600K we 
were able to make the assignment of the observed CO) lines 
without encountering serious difficulty, while at 800K an 
automated Loomis~Wood diagram technique was developed for the 
purpose of simplifying the identification. 

An absorption cell which could be operated at elevated 


temperatures was designed and constructed. The instrumentation 
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problem was mainly a question of stability in the optical 
path during high temperature operation. Prior experience 
indicated that the White cell configuration was totally in- 
adequate for high temperature operation and that a Pfund 
cell configuration might work very well. Even so, the multi- 
pass optical system was found difficult to Saeuate at high 
temperature - it worked within a limited temperature range. 
The mirror surface coated with a Rh film was adequate for 
the measurements at 600K and 800K. The cell transmission was 
found to deteriorate rapidly as the temperature was raised 
above 1000K. A single-pass optical configuration which 
required no mirror optics seemed the only arrangement with 
promise of successful operation at temperatures above 1000K. 
The measurements consisted of obtaining absorption data 


of both gases at two temperatures, 600K and 800K in the spectral 


region between 1600 to 2500 em, Analysis of the CO, data 


clearly indicated that a high accuracy of the measurement of 
the line positions is inadequate for extrapolating to the 
highest rotational transitions. The vibrational-rotational 
transitions of a linear polyatomic molecule (for example the 


Q branch) can be given by the simple expression: 


o(J) = Gt + BY (s+) — D'fa'(T'+1)}2 + HY (a'(s'+1)]3 


- G" = B™ o"(o"+1) + D"™[I"(I"+1)]2 - AM [I"(I"+1)] 3. 


It was found that the parameters, B, D, and H, determined from 


the measurement had a meaningful accuracy only within the 


measured transition range. They were by no means effective 


in predicting the transitions outside the range. For example, 


we found that those parameters determined from the data 
between P(50) to R(50) with high accuracy could predict P(60) 
to R(60) lines to poor accuracy. 

The intensities of the CO) lines observed in the data 
were not determined absolutely. In complex spectra such as 
those reported on here, very few observed spectral features 
are composed of single transitions. In dealing with the 
complex data we must determine the individual line intensities 
simultaneously for multiline components. The assignment of 
the exact positions for the component lines is of crucial 
importance to the effort. Because of this, the identification 
of all transitions observed in the data is of prime importance 
for. the data analysis. An extended effort is required for 
the assignment of transitions. 

The H20 data contrasted in many ways with the CO) data. 
The H2O molecule is an asymmetric top rotor; the rotational 
lines are distributed ina manner not so simple as a linear 
molecule and its vibrational frequencies are much larger 


than those of COQ). Even though the H2O spectrum contains 


far less lines than the CO) spectrum, its structure is far 
more complex than the latter. Nevertheless the analysis 

went well. The identifications of the transitions for H20 t 
and COQ2 and the intensities for H20 were made using the 


600K and 800K data. 
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It. GENERAL CONSIDERATIONS 

In all our observed spectra the 100% transmittance level 
was not constant but varied as a function of wavenumber. 

This variation was due to several factors such as, a variation 
in detector spectral sensitivity, the spectral transmission 
characteristics of the bandpass filter and probably of the 
interferometer as well. By far the most bothersome background 
level modulation was caused by the light interference between 
various optical surfaces which the infrared beam traversed. 
The molecular absorptance (for the case of H20) was determined 
with respect to the non-uniform 100% transmission level and, 
in spite of using various computation logics to reduce the 
effect of the background uncertainty, the uncertainty in the 
absorption was still greater than 5%. 

The absorption spectra, determined as indicated above, 
were then analyzed to obtain transition frequencies. Since 
the computer output is a spectrum consisting of discrete 
values in wavenumber at an interval corresponding to the 
reciprocal of twice the maximum path difference in the inter- 
ferogram, it is not smooth enough for analysis, i.e., connecting 
the spectral output points by straight lines results ina 
very jagged spectrum. It is necessary to determine extra 
spectral values between the original output points either by 
zero-extending the interferogram and Fourier transforming, 
or by convolving the original spectral output with a sinc 


function. Once this is accomplished it is easier te proceed 
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with the analysis. 


The line centers were determined either by using the 
calculated first and second derivatives of the smooth spectrum 
or by a technique described in Sec. IV. The uncertainty in 
the spectral line center was affected by the noise in the 
data as well as by the blending of neighboring lines. This 

é uncertainty is larger for the lines that are strongly saturated 
at the line centers, even for the case when these are well 
isolated from their neighbors. The overall uncertainty in 
the measured transition frequencies was close to the separation 
of the spectral points in the interpolated spectrum, i.e., 
about 0.0004 cm-l, 

The absorption contour of a molecular transition in a 
measured spectrum may be expressed, to a good approximation, 
by a convolution integral of the true absorption contour and 
the impulse response (commonly called the instrument function) 
of the instrument. For the spectrum obtained by using the 
technique of Fourier spectroscopy, it is more conveniently 
expressed in the interferogram domain (defined in the esticai i 
path difference scale x) rather than in the spectral domain k 
(defined in the wavenumber scale gg). The observed contour 


A'(a) and the true contour A(o) are related by 


-1i2n0x i2nox ; 
fA'(o) e do = {fA(a) e do} T(x), (1) 


where T(x) is a multiplicative function of finite extent 


limited to the maximum optical path difference X in the 


ee 
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; interferometer. The function T(x) has a shape which depends 


on the type of apodization applied to the spectral data. In 


our case, it is a triangular function tapering to zero at the 


maximum path difference X; 


Ix] 
T(x) = 1 - om 
xX 


if Ixl < 


and = 0 if Ix[| > x. (2) 


The true absorptance contour of a well-isolated molecular 


transition line is expressed by 


A(o) = 1 ~- exp [-k(o)], (3) 


where the function k(g) is the absorption coefficient defined 


by the transition strength S per single molecule in the form 


J(o0) do = SN. (4) 


In this formulation we assume that the absorption is caused 


by a uniform column of N total molecules per unit cross-sectional 


area, If thermal equilibrium exists along the absorption 


path, the strength S is given by the well-known formula 


hy 
8xv — 1 
gS: —- |RI2 (l-e kT) - e~E"/kT, 
3he g 


ey penne! 


where v is the transition frequency in Hz, R is the transition 


a 


moment, g is the statistical weight of the lower state, and 


E” is the lower state energy. The absorption coefficient 


k(oj takes on various shapes in accordance with the line 


profile. If the assumption is valid that collisions from 


neighboring molecules dominate the spectral line shape, we 


May express the absorption coefficient by the Lorentzian 


profile: 


Sa 
ko) > oy (5) 


n{ (0-09) 2 + a] 


where qa is the width and oo is the center frequency of the 


line. 


The observed spectrum for most cases contains lines 


which are to a varying degree overlapped with their neighbors. 


The absorptance contour for the overlapped lines is given by 


N 
A(o) = 1 - exp - [£ Knlo)), (6) 
n=] 


where the absorption consists of N lines, and kpy(o) is the 


absorption coefficient for the n-th line. The data which we 


deal with in the analysis contain all possible situations, 


i.e., from the case of well isolated lines to that of many 


overlapping lines. The objective of the spectral analysis 


is to determine the line parameters of each transition, the 


transition frequency oo, and for water the strength S and the 


width a. We describe below the technique used in extracting 4 


these parameters for overlapping lines. 


First, we separated the data into groups of lines which 


were sufficiently isolated and then proceeded to apply a 


curve fitting technique based on the least-squares method. 

The dapavdtion oF the lines into groups was necessary because 
the data contained lines of varying degree of overlapping. 

A criterion was set which established when the absorptance on 
either side of a line was sufficiently close to zero; the line 
was then considered to be isolated. With the line center 
frequency determined as described above, the remaining two 
variables for each component line (for the case of water), 

the strength and the width, were determined under the assumption 
that the line profile was Lorentzian. A spectral pattern 
‘theoretically constructed using the Lorentzian profile with 
assumed parameters was compared with the observed data. The 
square error between them was calculated and then minimized 

by adjusting the two variables, the strength and the width, 


for each component line. The sequence, which was the non- 


linear least-squares curve fitting process, was repeated to 
reduce the square error to an expected stationary value 

according to the noise level in the observed data. Once the : 
error figure converged to a stationary value, the integrated 


absorptance W, given by 


5 he Citi 


W= fA(o)do, (7) 


was calculated for each line contained in the group using 
the strength and the width thus established. The standard 


method of obtaining the integrated absorption by numerically 


integrating the observed absorptance for each line was not 
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applicable to our data reduction because the data contained 
considerable overlapping of lines. 
III. INSTRUMENTATION 

All the measurements were made using the Air Force 
Geophysics Laboratory (AFGL) high resolution Fourier transform 
spectrometer with the associated Pfund configuration hot 
cell. The maximum path difference the interferometer was 
driven to was 75 cm, resulting in an unapodized resolution 
of approximately 0.007 em7l, The AFGL step and hold inter- 
ferometer (built by Idealab) and hot cell system have been 
described in previous publicationsl.2+3, so will not be des- 
scribed here. However, we would like to indicate one modifi- 
cation in the operation of the system. 

In the past, the stepping of the interferometer mirror 
was controlled with dedicated electronics. With stable 
external conditions and the gains properly tuned, this hardware 
stepping system was adequate to step and hold 106 steps. 
However, in practice these conditions were found difficult to 
achieve and maintain. Consequently, a PDP-8E minicomputer has 
recently been used to increase the stepping reliability of 
the interferometer. This mini-computer, used previously only 
for handling data, is now used to monitor the laser reference 
signal and take corrective action for the most commonly 
occurring stepping errors. In addition to relaxing the con- 
ditions under which the interferometer can be successfully 


operated, this computerized stepping control system provides 
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very useful diagnostic information on the operation of the 
stepping system. Before the computerized control system 
was added it was almost impossible to determine anything 
about the nature of stepping errors that occurred infrequently. 

For the actual measurements the high temperature absorp- 
tion cell, inside the oven, was placed against the entrance 
window of the interferometer vacuum chamber. Since we were 
unable to properly couple the hot cell to the interferometer 
chamber there remained about a 0.25 cm optical path through 
ambient air. The remainder of the optical path from the source 
to the interferometer was maintained below 1 torr; except of 
course when the gases were introduced into the absorption 
cell when the pressure then was either 3 or 6 torr. The 
pressure in the absorption cell was measured by a Baratron 
pressure guage. 
IV. CO MOLECULE 

The importance of CQ) as an atmospheric molecule justifies 
measurements and calculation at high precision of band para- 
meters for a great number of bands; it is indeed a significant 
contribution to both theoretical and experimental studies 
of our atmosphere. Very accurate parameters on the (00001-00011)* 


band of 12¢16q, have been obtained by Pine and Guelachvili’. 


*For CO2 there are three fundamental modes of vibration, vy, v2, 
and v3. Associated with the bending mode, v2, is the angular 
momentum ¢. The notation used is v)votv3r where r is the ranking 
index for a Fermi resonating group. 


Their measurement was accomplished by combining two spectra, 


a room temperature Fourier transform spectrum and a high 
temperature tunable laser spectrum. There have been numerous 
Measurements of other bands of C02, but they were made either 
with lower resolution spectrometers, near room temperature 


where fewer rotational lines were excited, resulting in a 


pba atts 2a Lae meena eet Si 


less accurate determination of band parameters, or electrical 
discharge spectra where high vibrational levels were obtained 
but with the normal distribution of rotational level J. In 
this section we present measurements of CO) made with the 
above described system. 

The energy (in em71) of a linear molecule, such as COo, 
in a given vibration-rotation state can be represented by 

E(v,J) = G, + By(J(J+1)] ~ D,(a(s+1)]? + Hy[a(g+1)}3 
where v is the quantum number associated with the vibrational 
energy, J is the quantum number associated with the total 
angular momentum of the molecule, and G, B, D are the spectro- 
scopic constants. The infrared spectrum corresponds to 
transitions between different vibration-rotation states. 

The advantages of making the measurements at elevated 


temperatures are twofold. First, at higher temperatures more 


“rotational lines are excited, making it possible to obtain 


t 
{ 
i 
higher precision in the determination of the spectroscopic t 
constants G, B, D, and H. Secondly, the higher temperature 


results in higher vibrational energy levels being excited, 


making more of the so-called “hot bands" measurable. " 
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We used a PDP-8e mini-computer to control the inter- 
ferometer and to sample the interferogram. The data were 
taken via magnetic tape to be analyzed at AFGL and/or at the 
University of Massachusetts. The data were first Fourier 
transformed, the various lines identified, and finally a 
least-squares fit technique used to generate new values for 
the band parameters. Each band was fit independently of 
the other bands. 

Spectra obtained using high resolution Fourier trans- 
form spectrometers have an excellent frequency stability 
over a wide spectral range. There is, however, a small 
systematic frequency shift introduced into the spectra due 
to finite detector size. In principle, this correction can 
be calculated from the geometry of the interferometer apertures 
and detectors, but in practice it is usually easier to use 
an internal frequency standard such as lines of CO which 
were present in our experimental spectra. Guelachvili°® has 
measured the positions of these CO lines with an accuracy of 
0.00008 cm-1. We used these line positions as an internal 
frequency standard. After this systematic error was removed 
from our spectra, the standard deviation of the fit between 
our data and Guelachvili's, for well isolated lines, was 
about 0.0004 cm-l, 

The identification of the transitions was made in one 


of three ways, depending on how well the band had previously 


been measured. For well-known bands the identification ra 


; 
I 
i 


22 


at 


procedure was readily automated. The position and strength 


of each line was calculated, and then used to locate the line 
in the experimental spectrum with approximately the right 
strength that was closest to the calculated position. 

For other bands, the positions of lines for low J have 
been well measured, but not for high J. In this case, the 
line positions were again calculated and the experimental 
lines identified, starting at low J and moving up to higher 
J, until there were not any lines within about 0.02 cm! of 
the calculated line. The band was then refit and new constants 
obtained. The process was iterated until no further extension 
to higher J was possible. 

In some cases the band was so poorly known even at 
low J that the identification could not be made by taking 
the experimental line closest to the calculated line position. 
Another problem arose when certain combinations of merged 
lines made the extension from low J to higher J impossible. 

In these two cases, a Loomis-Wood diagram® approach has been 
extremely helpful in picking out the lines that belong to 

one band in the presence of lines belonging to other bands. 
When the differences between a calculated line and all observed 
lines that are nearby are displayed graphically, it is rela-~ 
tively easy to see the pattern created from a set of lines 
which have spectroscopic constants similar to those used to 


calculate the line positions. 


Iv.a. 12cl6o, (626)*, 2180 cm™1 to 2400 cm™!, 800K 


| 

The data presented in this section were obtained at 800K | 

with a resolution of 0.007 cm~1. once the identification | 

of these measured lines had been made each band was fit ' 

discarding only the most severely merged lines. Lines for 
which the difference between the observed and calculated 


position was less than about 0.005 cm] were kept in the fit. 
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All lines used in the fit were weighted equally. The bands 


that were fit are listed in Table I, along with the standard 


| 

\ 

‘ deviation and the range of J values used. The resulting 
values obtained for G'-G", B', D', H', B"™, D" and H" are 
tabulated in Table II. 

Several constraints were used on the fits. The spectro- 

scopic constants H' and H” were constrained to zero for some 

bands. For bands where zg is not equal to zero, r-type doubling 

| results in two sets of energy levels denoted by c and d. 

| The spectroscopic constants obtained for the c and the d 

| levels were considered separately for bands where g=1 and 2=2. 
For the bands where G" is less than 3000 cm7! restrictions 


were placed on the G's and B's. For g=l the G's, and 


for e£=2 both the G's and the B's of the c levels were 


The difference between the transitions calculated using 


3 
| constrained to be equal to those of the d levels. 
the band parameters determined by Guelachvili and those in 


the 1978 AFGL line compilation® for the P branch of the 


* Designation used on AFGL tapes’ 


' \ 
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Table I. The 12¢l6o, (626) Bands Observed at 800K 


Transition Band Origin Range of Standard 


(cm7}) Measurement __ Deviation (cm71) 


00011 00001 2349 .143 P120-R118 0.0009 
01111¢C 01101C 2336 .633 P107-R109 0.0007 
01111D 01101D 2336 .633 P108-R104 0.0006 
10012 10002 2327 .433 P102-R102 0.0008 
02211C 02201C 2324.141 P98-R104 0.0009 
02211D 02201D 2324.141 P105-R107 0.0009 
10011 10001 2326.598 P104-R104 0.0009 
03311 03301 2311 .668 P90-RI2 0.0009 
11111¢ 11101C 2313.772 P89-R87 0.0009 
11111D 11101D 2313.772 P88-R88 0.0012 
20013 20003 2305 .257 P80-R76 0.0013 
12212C 12202C 2302 .958 P68-R62 0.0026 
12212D 12202D 2302 .958 P67-R67 0.0026 
04411 04401 2299 .214 P82-R83 0.0010 
20012 20002 2306 .690 P86-R56 0.0013 
12211C 12201C 2301 .053 P86-R88 0.0010 
12211D 12201D 2301 .053 P83-R83 0.0010 
00021 00011 2324 .183 P83-R87 0.0007 
20011 20001 2302 .523 P82-R78 0.0011 
01121C 01111C 2311.703 P90-R80 0.0013 
01121D 01111D 2311.700 P83-R79 0.0014 
10022 10012 2302.369 P69-R65 0.0018 
02221C 02211C 2299 .240 P75-R93 0.0013 
02221D 02211D 2299 .238 P78-R80 0.0012 
13311 13301 2288 .390 P65-R63 0.0014 
13312 13302 2290 .680 P73-R66 0.0013 
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02201C - 02211C band of 42c16o, is shown in Figure 1. 
Guelachvili made extremely accurate room temperature measure- 
ments (0.0001 cm!) of the positions of transitions up to 
about J = 46 of this band, from which he calculated spectro- 
scopic constants?. For J < 50 the agreement with our data 
is excellent, but for high J the difference is large. 

The large values of the standard deviations given in 
Table I, considerably greater than the 0.0004 cm7! for well 
isolated lines mentioned earlier, are caused by merging from 
overlapping bands. Most observed lines are merged to some 
degree. If the overlapping lines belong to bands that are 
not correlated with the band being fit, the effect of the 
merged lines will be to produce an increase in the random 
scatter of the fit. It is very rare to have both the position 
and the spacing between the lines of the band being fit and 
the lines in a merging band close enough that significant 
correlation occurs. The strongest correlation between bands 
in the observed spectrum appears to be between c and d bands. 
For example, the 01101C-01111C and the 01101D-01111D bands 
have spectroscopic constants just different enough that even 
though they have the same band origin, R(105) of the c band 
is merged with R(106) of the d band. At this point in the 
bands, the difference between the line spacing is only 0.005 
em7l, yet R(103) is not merged with R(104), nor is R(107) 


merged with R(108). The smallness of the correlation between 
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-ands supports the assumption that fitting partly merged 
lines does not introduce systematic errors. 
To further insure that the identification and fitting 

procedure was valid, a synthetic spectrum was calculated and 

3 compared to the observed spectrum. The spectroscopic constants 
listed in Table II and the band intensities given in reference 
7 were used to create the synthetic spectrum. Bands that 
resulted in higher absorption in the synthetic than in the 


observed spectrum were either refit or thrown out entirely. 


13,16, 13.16.18 


IV.b. 2 and ““c™°0""0; 2140 cm”? to 2340 om7?, 600K 


Spectroscopic constants for L3cl6o, and 13cl6ol8o that 
predict the position of spectral lines for J values greater 
than 100 were obtained from absorption data of an isotopically 
enriched sample of C02. The gas sample, consisting of 888% 
13¢l6q,, 11% 13¢cl6ol85 ana 1% various other isotopes, was 
heated to a temperature of 600K. By making measurements on 
several isotones a more accurate determination of the electric 
potential for the CO) molecule can be made since the electric 
potential function is identical for different isotopes of 
the same molecule. 

Each spectrum was computed by averaging three transformed 
interferograms, each of which had a measurement time of 
fifteen hours. Spectra were taken at both 300K and 600K, 
but only results of analysis of the 600K spectrum containing 
high J lines are presented. The 300K spectrum was taken in 


order to assist in the identification of 13c160180 bands in 
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the 600K spectrum. 

To eliminate the effects of a sloping background, a 
spectrum was taken of the empty absorption cell and the 
desired spectrum was ratioed against this spectrum. To 
minimize the effects of noise in the empty-cell spectrum a 
13 point running average was first used to smooth the spectrum 
before the ratioing was performed. 

The spectrum was then analyzed to obtain the position, 
asymmetry, width, and strength of each line. The density of 
the spectral points was increased by a factor of sixteen 
through interpolation. The line positions were determined 
by taking the wave number position of the point closest to 
the local minimum as the line center location. A measure of 
line asymmetry was given by using an alternative method for 
finding the line positions and comparing the results of the 
two methods. This second method used the center of a chord 
drawn across the absorption line a small distance up from 
the minimum to determine the line center. For symmetric lines, 
the center of this chord coincides with the line position 
found from the local minimum method. The line width was 
determined as being proportional to the length of this chord. 
An estimate of line strength was calculated by integrating 
the area between the spectral line and the background. The 
asymmetry, width, and line strength were used to determine 
the amount of line merging present and hence the quality of 


each line. 
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As a check on the quality of the experimental data a 
comparison of the line positions of the 00001 to 00011 band 
of 12¢l6o,, observed from P(72) to R(78) was made with 
the line positions computed using the spectroscopic constants 
given by Pine and Guelachvili4. ‘The results of this comparison 
showed a systematic difference of 0.0001 em-l and an rms error 
of 0.0008 cm-l. 

For low J lines, identification was accomplished by 
assigning the line in the experimental spectrum whose position 
was closest to the position calculated using the AFGL (1980) 
spectroscopic constants!l, the lines were identified 
starting oe low J and moving up to higher J until there 
weren't any lines close (within about 0.02 cm71) to the calcu- 
lated values. The band was then refit, new constants obtained, 
and the process iterated until further extension to higher J 
was impossible. For the v3 fundamental band of 13¢16o185 
it was not possible to make conclusive identification of the 
lines in the 600K spectrum, without first using the AFGL (1980) 
constants to identify the lines in the 300K spectrum. The 
band was first fit using the 300K data. The resulting con- 
stants were then used to make the identification in the 600K 
spectrum. 

After the lines were identified, a fit was then performed 
on each band. All lines within 0.0025 cm-! of the calculated 
positions were kept in the fit. In order to insure that the 


identification process was working properly, lines that were 
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not included in the fit were checked for asymmetry and abnor- 
mally large width or strength to verify that their exclusion 
from the fit was justified. All lines used in the fit were 
weighted equally. The bands that were fit are given in Table 
III along with the range of J values, the total number of 
lines, and the standard deviation of each fit. The resulting 
values obtained for G'-G", B', D', B", and D" are tabulated 
in Table IV. All spectroscopic constants were allowed to 
float freely in the fit. No attempt was made to combine these 
constants in such a manner as to yield a single consistent 
potential function for the CO) molecule. 

For the 00001 to 00011 and the 01101 to O1111 bands of 
13¢16o, the low J lines were badly saturated. Guelachvili 
had previously measured these lines to a high accuracy?, 
hence for low values of J the line positions of Guelachvili 
were used in the fit. Guelachvili's line positions were used 
for the P(52) to R(50) lines of the 00001 - 00011 band, from 
P(35) to R(25) for the 01101C - 01111C band, and from P(40) 
to R(40) for the 01101D-01111D band. 

Comparisons of the line position measurements determined 
in this work to those predicted by the extrapolation using 
Guelachvili's constants? and those used in the 1980 AFGL 
line compilationl! are given in Figures 2 to 4. 

Several interesting features are illustrated by these 


figures. Although Guelachvili's line positions were measured 


with considerably more accuracy than ours, extrapolation to 
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higher J (higher temperature) is unsatisfactory. It appears 


to be a general principle that spectroscopic constants obtained 


using low J lines do not successfully predict the position 
of high J lines. By comparing the difference in the quality 
of the AFGL (1980) fit of the 01101C-01111C band to the D 
band it can be seen that the AFGL (1980) line compilation is 
quite good for some of these very weak atmospheric bands 
while not for others. The spectroscopic constants for the 
eight bands reported in this section predict the position of 
spectral lines to an estimated absolute accuracy of 0.001 
cm-l for J values up to 100. This represents a significant 
improvement in the accuracy with which high J line positions 
of the isotopes 13¢16o, and 13¢160186 can be computed. This 
improvement illustrates the advantage of heating the sample 
to determine spectroscopic constants that are valid for 
large J values. For several of these bands a positional 


accuracy improvement was realized for low J lines as well. 


ree Os and '3¢!®q18q, 2140 cm) to 2340 cm 2, 800K 


IV.c. 
In this section we present results of measurements of 
the same isotopic species as the previous section but at a 
temperature of 800K. 
A total of four 800K spectra were taken, two with a 
CO2 gas pressure of 3 torr and two with 6 torr. In addition 
to these four spectra, several empty hot cell spectra were 


also recorded to be used as a reference. The measurement 


time for each spectrum was approximately fifteen hours. 
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After the two 3 torr spectra and the two 6 torr spectra were 
each separately co-added they were ratioed by the empty cell 
reference spectrum to reduce background effects. To minimize 
noise in the empty cell spectrum a thirteen point running 
average was first used to smooth the spectrum before the 
ratioing was performed. The 3 torr and the 6 torr spectra 
were then analyzed as described in Sec. IV.b. to obtain the 
position, asymmetry, width, and strength of each line. 

The asymmetry, width, and line strength were used to 
determine the amount of line merging present and hence the 
quality of each line. Only severely merged lines were excluded 
from the least squares procedure. Lines which were slightly 
merged were still used but with a reduced weighting. The 
weight associated with each line, the reciprocal of the expected 
uncertainty squared, was estimated from the line asymmetry, 
abnormal line width, and inconsistencies of line position as 
computed from other lines in the band. Before the final fit 
was realized the 3 torr and 6 torr data were combined to 
form a single set of data. 

Since there was a small amount of 12¢l6o, in the high 
temperature cell it was possible to make a check on the 
quality of the experimental data by comparing the line positions 
observed from P(66) to R(84) of the 00001 to 00011 band with 
the line positions computed using the spectroscopic constants 
given by Pine and Guelachvili4. The results of this comparison 


showed evidence for the existence of a systematic error, 
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particularly for the weak lines at the high end of the R 
branch where the systematic difference reached a maximum of 
0.0006 cm-!. This error probably resulted from incomplete 
phase correction. The average difference was less than 
10-4 cm-l with a rms difference of 0.0003 cm-l. 

A weighted least squares fit was then performed on each 
band. An example of the quality of the fit and the data 
used is given in Table V for the P branch of 13¢160, for the 
00001 to 00011 band. This table lists the observed line 
position, line identification, difference between observed 
and calculated position, and the expected uncertainty for each 
line. 

Since the low J lines for the 00001 to 00011 band of 
13¢16o, were badly saturated for both the 3 torr and 6 torr 
spectra, lines from P(40} to R(40) of this band were excluded 
from the fit. 

The bands that were fit are given in Table VI, along 
with the range of J values, the total number of lines used, 
and the standard deviation of each fit. The resulting values 
obtained for G'-G", B', D', H', B", D", and H" are tabulated 
in Table VII. All spectroscopic constants were first allowed 
to float freely in the fit. The resulting values of H' and 
H" were smaller than the uncertainty in these parameters for 


some bands, therefore H' and H" were constrained to zero. 


These bands are designated in Table VII by letting H' = 


“chk A sched chtmeasaal 


Table V. The P-Branch of the 00001 to 00011 Band 
of 13¢1604 (636) 


Observed 
Minus Expected 
Assignment Observed Calculated Uncertainty 
P(42) 2245 .6445 0.0001 0.0010 
ae 2243 .5859 0.0002 0.0008 
P(46 2241 .5038 -0.0001 0.0008 
te 2239.4004 0.0013 0.0017 
P(50 2237 .2715 0.0004 0.0005 
P(52) 2235.1204 0.0002 0.0009 
P(54) 2232 .9463 0.0001 0.0005 
bles 2230.7495 0.0002 0.0005 
} P(58 2228 .5298 0.0004 0.0005 
+ P(60) 2226 .2873 0.0006 0.0008 
, P(62) 2224 .0212 0.0002 0.0004 
a P(64) 2221.7328 0.0003 0.0005 
; P(66) 2219.4215 0.0003 0.0004 
P(68 2217 .0873 0.0001 0.0004 
P(70 2214.7288 -0.0015 0.0013 
P(72) 2212.3504 -0.0004 0.0005 
P(74 2209 .9480 ~0.0005 0.0007 
P(76 2207 .5236 -0.0000 0.0004 { 
P(78) 2205 .0754 -0.0007 0.0005 4 
P(80 2202 .6041 -0.0018 0.0012 F 
P(82 2200.1125 ~0.0007 0.0005 
P(84 2197 .5983 0.0003 0.0007 
P(86) 2195.0597 -0.0006 0.0005 
P(88 2192 .5003 0.0002 0.0006 
P(90 2189.9165 ~0.0009 0.0011 
P(92) 2187 .3105 -0.0019 0.0012 
P(94 2184 .6860 0.0010 0.0010 
P(96 2182 .0354 0.0002 0.0004 
P(98) 2179.3648 0.0016 0.0013 
P(100) 2176 .6696 0.0007 0.0007 
P(102) 2173.9520 -0.0003 0.0004 
P(104) 2171.2127 -0.0008 0.0011 
ie 2168.4526 0.0000 0.0003 
P(108 2165 .6698 0.0003 0.0004 
P a 2160 .0387 0.0016 0.0013 
P(114 2157.1878 -0.0000 0.0006 
aie} 2154.3160 -0.0906 0.0011 
P(118 2151.4217 ~9.0016 0.0023 
ioe 2148 .5085 0.0003 0.0011 
P(122 2145 .5735 0.0024 0.0024 
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For the bands where g-type doubling is present the c 
and d levels were fit independently. Since the band origins 
resulting from fitting the c levels should be identical to 
the one resulting from the d levels, a check on the precision 
of the data can be obtained by comparing the difference 
between the two. The difference between the band origins for 
the 01101 to 01111 transitions was 0.0003 cm-l, and for 
the 02201 to 02211 transitions was 0.0001 cm-l. 

Comparisons of the line position measurements determined 
in this work to those predicted by the constants used in the 
1980 AFGL line compilation and Guelachvili'ts constants for 
those bands which he measured? are given in Figs. 5 through 
12. Several interesting features are illustrated by these 
figures. Although Guelachvili's line positions were measured 
with considerably more accuracy than ours, extrapolation to 
higher J (higher temperature) was unsatisfactory, as shown 
in Figs. 5 through 7. However, interpolating from high J to 
low J seems to work quite well as indicated by the comparison 
shown in Fig. 5 between the fit of the 00001 - 00011 band 
of ed ooh excluding lines from P(40) to R(40), and the 
measurements of Guelachvili. The difference in the band 
center as compared to Guelachvili's for this band was only 
0.0002 cm~1. Table VIII is a listing of the lines used in 
the fit. 


By comparing the difference in the quality of the AFGL 


(1980) fit of the 01101C-01111C band to the D band (Figs. 6 
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Figure 5. Comparison of measured line positions with those 
computed using Guelachvili's and the AFGL (1980) 
constants for the 00001 to 00011 band of 13cl605. 
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Figure 7. Comparison for the 01101D to 01111D band of l3cl€o,. 
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Figure 9. Comparison for the 02201C to 02211C band of 13¢160,. 
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Pigure 10. Comparison for the 02201D to 02211D band of 13cl6o,. 
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Figure 12. 
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Comparison for the 10001 to 10011 band of 13cl6o05. 
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and 7) it can be seen that the constants used for the AFGL 
(1980) line compilation are quite good for some of these 
very weak atmospheric bands but not so good for others. 


The improvement in the spectroscopic constants obtained 


at 800 K over those of the previous section for 600 K are 
not particularly significant except for the 02201C to 02211C 
band of 13¢l6o, where the P branch shows a marked difference 
(Fig. 9). Using the expansion in powers of J(J+l) (Eq. 1) 

to extrapolate from a lower to a higher J is unsatisfactory, 
particularly for this band as indicated by the poor positional 
accuracy of high J lines predicted using both the AFGL (1980) 
constants and the 600 K constants. Also, there is a tendency 
for oscillations in the fit. The H values obtained for this 
band are unusually large and negative. Large negative H's 
were also reportedl2 (table II of Sec. IV.a.) for these 

same two vibrational levels 02201C and 02211C for the main 
isotope 12¢16o,. 


The spectroscopic constants for the eight bands reported iy 
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in this section predict the position of spectral lines to an d 
estimated absolute accuracy of better than 0.001 em-1 including 
J values greater than 100. This represents a significant 
improvement in the accuracy with which high J line positions 
of the isotopes 13¢16o, and 13¢l6o186 can be computed. 
V. «B90 (161); 1600 cm! to 2001 cm™1, 800K 

In this section we present results of measurements of 


H20 using the instrumentation and the analysis technique 


rae 
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described above. 

A Lorentzian profile was assumed as the line shape in 
our analysis - which is probably improper in principle. 
However, the H20 data were always taken at very low pressure 
in order to reduce to a minimum the line overlapping problem. 
The expected collision width was much smaller than the Doppler 
width under such a condition. The Lorentz line shape was 
adopted because of the following reasons. For weak lines 
the question of line width is irrelevant because the inte- 
grated absorptance W and the strength maintain a linear 
relationship independent of width. For strongly saturated 
lines, the observable absorptance contour corresponded to 
the wing region of the line which was Lorentzian. The central 
peak region of such lines was completely saturated. The 
integrated absorption W was controlled by the wing section 
and thus by the Lorentzian profile. The line profile in the 
center peak region produced only a small effect in the deter- 
mination of the integrated absorptance. Thus the Lorentzian 
profile assumption was a proper practical choice for obtaining 
the integrated absorptance even though it was theoretically 
improper. 

Once the integrated absorptance was determined, it was 
used to derive a value for the final strength. For the 
Lorentzian line which was strongly saturated at the center, 
the strong interdependence existing between the strength S 


and the width aq made their separate determination very 
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impractical. With the integrated absorptance and the absorp- 
tance contour fixed, a quantity /Sa was the only meaningful 
parameter which could be determined for such a line, even 
when the noise level in the data was exceptionally small. We 
adopted the theoretically accepted value for q (Ref. 13) in 
deriving the strength from the integrated absorptance. The 
width was extrapolated at 800K from the values calculated 
for room temperature using a simple impact theory assumption; 
T) 1/2 
a (T) = a (T)) (- ) Z 
T 

The results obtained for water vapor are summarized in 
Table IX. The data listed in the table are the line frequency 
in cm-l, the intensity S in (cm-l/molecule cm-2) at 800K, the 
observed integrated absorptance for an observation condition 
PHO = 6.0 torr, T = 800K and zg = 350 cm, the lower state 
energy in cml, the transition identification (J', Ka', 
Ke'), (J", Ka", Ke"), (v'y, v'o- v'3) and (v"y, v"2- v"3), 
and the difference between the observed transition frequency 
and the value listed in the latest AFGL line listingll, 
The data are arranged in ascending order of observed line 


position. The data listed in Table IX contain those lines 


ove ET Sie eon 


newly identified in the present study, which are the high J 
lines of the (v2-0) and the (2v2-v2) transitions. We were 
able to follow the v2 transitions up to the (241 ,24-239,23)- 


(249 ,24-231,23) doublet, and the (2v2-v2) transition up to 


the (20), 2071909,19)-(209, 20-191,19) doublet. The highest 


excitation energy observed at 800K in our data exceeded 5000 
cm-!, Table X summarizes those lines which were either 
newly assigned or quite different from the latest AFGL line 
listing. 

VI. Global Constants for Carbon Dioxide. 

One of the principal aims of the research described in 
this report has been to provide data from which reliable 
parameters, necessary for high resolution atmospheric trans- 
mission/emission calculations, could be obtained. In this 
section we describe the application of the AFGL FTS-hot cell 
measurements on CO) to the determination of "global constants". 
This approach is to be distinguished from the discussion in 
Sec. IV, where vibration-rotation bands were least-squares 
fit independently to obtain spectroscopic constants that 
would best reproduce that particular band with the minimum 
root-mean-square deviation, regardless of the fact that the 
vibrational levels involved in the transition are themselves 
part of an ensemble of many vibrational levels connected by 
allowable electric dipole transitions. In the approach for 
global constants, all connected vibrational states that have 
been observed at high resolution (in the laboratory) are 
taken in concert. The approach is analogous to forming a 
“tree” and bootstrapping to successive vibrational levels. 
For each isotopic species we thus require the constants to 


form a self-consistent set. In general the deviations between 


observed and calculated positions of transitions will now be 
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Slightly larger than the best fit of individual bands, 
depending on the weakest link in the path, either the accuracy 
of the measurement of a band, or more usually, the limit of 
the highest value of rotational quantum number observed. 
Nevertheless, a great consistency has been realized by this 
method for generating the line positions of the approximately 
560 vibration-~-rotation bands of importance to atmospheric 
optics problems in the IR, covering a spectral range of 400 
to 10,000 cm-l. the method has been discussed somewhat in 
Ref. 10, but several points will be elaborated upon here. 

It should be remarked that this method does not go to the 
extent of combining the various isotopic measurements com- 
pletely to form a true valence bond force potential, (see 

for example ref. 14) but treats them quasi-independently, 
only making extrapolations for the other isotopes where 
measurements have not been made. Common to both approaches, 
as has been discussed in Sec. IV, is the difficulty in making 
accurate predictions of rotational levels much beyond the 
highest measurements. It is indeed this aspect that has 

made the results of the measurements of the high temperature 
CO2 at AFGL such a step forward in the determination of good 
parameters since the calculation of intermediate rotational 
levels is reliable and hence the extension of J-values that 
has been described here has greatly extended the reliability 
of line parameter compilations. 


We first review some of the considerations involved in 
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pian Scataccgacs: 


a description of the energy level and rotational constants 

for CO2. For a linear triatomic molecule like C09, the 
vibrational states are characterized by three quantum numbers, 
which are zero or positive integers, of pure vibration, vj, v2, 
and v3 and a fourth number 2 = v9, v9~-2,..- -v2 + 2, -v2 which 
represents the contribution of the bending mode to angular 
rotation. In order to calculate the purely vibrational part 
of the energy, it is necessary first to compute an unperturbed 


energy, 


unp . 2 
Gy P= ) wiVq + Y *i5%i93 + 922 2° + e: YijkViYj"k + 
1 13 ijk 


(8) 
) Yiee Vi a 
1 
and then to incorporate the effects of resonance perturbation 
by combining all close-lying levels with common ¢g and common 
symmetry in matrices whose diagonal elements are G,"P, whose 
off-diagonal elements are functions of additional molecular 
constants and the four quantum numbers, and whose eigenvalues 
are the vibrational energy Gy. 

This perturbation calculation results in the "mixing" 
of states whose G,""P are particularly close, so that the 
final description of the level by the original four quantum 
numbers is a poor one. Accordingly, it is useful to add a f 
fifth index, the rank symbol r to label all such mixed states if 
in order of decreasing energy. We have adopted this procedure 


for C02, where the mixed states are (vj) v2 £2 v3), (v1 + 1, 


‘ 
; 
( 
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V2 - 2, £27 V3)e + © » etc. The highest value of v, and the 
lowest value of v2 in each set are retained in the vibrational 
identification for all levels of the set. 

There is an accidental (Fermi) resonance between the 
vibrational states v] and v2 = 2v,. The wave functions 
become a mixture of each of the two states contributing to 
each perturbed level. 

We describe the energy levels by Viv2zv3r, where the 
ranking index, r, is unity for the highest vibrational level 
of a Fermi resonating group. The ranking index can assume 
the values l, 2, .. vy + 1. The quantum numbers v2 and ¢ are 
always equal in the AFGL notation. 

It will be noted in Eq. (8) that the vibrational energy 
depends on e2. When xz # 0 there are two levels for each 
J > & and this degeneracy is removed by rotation. The 
splitting ("xg-type doubling") results in two sets of 
levels, designated c and d, with different effective rotational 
constants. When g = 1 the splitting is most important, and 
Bo # Bar Do # Dg, etc.; when 2 = 2, Be = Bg but Do # Dg, etc.; 
when 2 = 3, the constants other than H are the same. These 
results are rigorous and arise from the effect of removing 
the degeneracies. For the pi-states (2 = 1) one performs a 
single contact transformation and the off-diagonal terms 
have the effect of a term in J(J+1), i.e., change the effective 


B-rotational constant. For delta-states (2 = 2), the coupling 


between diagonal numbers is twice removed, necessitating a 


a ey arenes a ea 


iil i is 


double contact transformation which has the effect of J2(J+1)2. 
This first affects the fourth order or D-centrifugal stretching 
constant. Resonances occasionally cause larger deviations, 
such as might be seen in some of the Hy constants. 

When the linear molecule has a center of symmetry, as 
in CO, with 165 at both ends (but not when one oxygen is 
isotopically different), the paired atoms with zero nuclear 
spin cause zero statistical weight for rotational levels of 
a given parity. Thus, only even-J levels exist for the ground 
vibrational level and for all other levels with ¢ = 0 and v3 
even ae symmetry"); for levels with zg = 0 and v3 odd (roe 
only J odd exists; when ¢ > 0, the c- and d- sublevels have 
different symmetry, so that for g = 1, v3 = 0 (My) the J = 
odd levels are c and the J = even levels are d, etc. In the 
AFGL compilation, symbols c or d are appended to the rotational 
quantum number of the lower state only when required, that 
is for g > 1 in the molecules without the center of symmetry. 
For example, R27C means 28c <« 27c; Q27C means 27d <« 27c. 
Figure 13 illustrates this scheme for point group Cay for 
the first few lines of a f-l transition.* 

The line frequencies are determined from the energy 
states by taking the differences corresponding to all allowed 


transitions. These depend on the familiar selection rules 


*T€ should be noted that the c and d labels correspond respec- 
tively to the new standardized e and f labelling of parity 
doublet levels (see Ref. 15). 
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Figure 13. 


P-branch Q-branch R-branch J 


+ 


+ 


+ 


Point Group C,, (for ‘OC’0) 


Designation of c- and d- levels for II-ll transition 
for COz without center of symmetry. 


for the linear molecule: 
when 4 £= 0,47 = +1, c*#c, d«d, 
when A 2 = 1,40 = +1, c= ¢, dd, 
and also AJ = 0, c « d. 

The line positions may thus be conveniently represented 
for computational purposes as a given series in m, where m = 
J* + 1 for the R-branch (J" + J"+1), m = -J" for the P-branch 
(J" » J"-1), and m = J" for the Q-branch (J" + J"). A different 
equation is needed for the Q~branch of a given transition 
than for the P and R branches, because of the differences 
in the c and d constants. The general equation is 


24 cm3 + am4 + em> + fm®, 


o = dg + am + bm 
where oo is the band origin, a = B'y + B"y, b = B'y = B"y + 
D"y - D'y, c = -(2D'y + 2D"y - H'y - Hy), d = Dy - 
D'y + 3H'y - 3H"y, @ = 3H'y + 3H"y, and f = H'y = Hy. 
For the Q-branch, the line positions are given by 
0 = gg + IN(I" + 1)(BYY — BMY)+I"2(9" + 1)2(D", = DIV) +... 
When expressed in this form it becomes clear from com- 
parisons between new observations and calculations which 
terms are contributing. to the errors. From the figures 
presented of the new data on CO, it can be seen that the 
even terms are most noticeably in error. This demonstrates 
the sensitivity of the experiment to the difference of con- 
stants such as the upper and lower B-values, or inverse 
moment of inertia, which are nearly the same for each state, 


The absolute values of the constants come more directly into 
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the odd terms in the running index m. For an adequate 
determination, it is thus necessary to use the global approach. 
The method employed in obtaining global constants has 
utilized the observations made at AFGL at high temperature 
from 1600 to 2500 cm-l in conjunction with high resolution 
measurements in several other key regions for CO). In the 
bending mode fundamental region (around 15ym) we have used 
the FTS measurements of Paso et all®, tn the CO laser region 
(~» 10um) we have incorporated the results of Freed et all17, 
Dupre-Maquaire and Pinson!8, and Siemsen et al!9. At shorter 
wavelength (» 2ym) we have used the high resolution measure- 
ments of Maillard et al29, In the same region that our 
interferometer has scanned, additional measurements at high 
vibration have been supplied to us from electric discharge 
measurements made by Bailly et al21. the band constants 
obtained from this effort will appear in a forthcoming publi- 


cation. 
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Table IX. HO Line Data Between 1600 cml and 2000 cm71 


The listings are: 
the transition frequency in cm-l 
the absorption strength in cm-l/molecule cm~2 at 800K 


the integrated absorptance in em~! for the H 0 column 
density 2.17 x 1020 molecules/cm2 at 800K 


the lower state energy in cm-l 


the transition assignment in (J', Ka', Ko"), 
(J", Ka"s Ko"), (vi'ys v2's V3"), and (vy", V2", V3") 


an internal code 
the isotope code 


the molecule identification code 


ee te ee Ee nd toe a FER EE tt Sart Beemer ate 


the integrated absorptance in em7! for the H5,0 column 
density 2.17 x 102 molecules/cm2 at 800K, and 


the wavenumber difference between the measured and the 
one listed in the AFGL (1978) listing. 
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Table X. 


Newly identified lines of the (010-000) and (020-010) bands. 


Transition Position 


(010-000) band 


211,21 - 209,20 


1949.9147 cm! 


219,21 ~ 291,20 


221,22 - 219,21 


19634065 4 


229,22 - 211,21 


231,23 - 229,22 


1976.7223 


239,23 - 221,22 


24) 24 - 239,23 


1989.9184 


249,24 - 231,23 


(020-010) band 


121,12 119,11 1780.9916 


129,12 111,11 1780 .9407 


13],12 - 129,1 1849.5795 


131,13 - 129,12 17970532 


139,13 121,12 1797 .0276 


141,13 139112 1868 .0298 


139,13 1812.8915 
blended 
131,13 1812.8915 


141,14 


149,14 


153,13 - 142,12 1931.9649 


152,13 - 143,12 1928 .6794 7 


wv 252,14 - 141,13 1888.0332 


* $15) 14 - 142,23 1887.7405 


Table X. 


(Continued) 


Transition 


151,15 


159,15 


163,14 
169,15 
163,15 
161,16 
169,16 
171,17 
179,17 
181,18 
189,18 
191,19 
199,19 
201,20 


209,20 


149,14 
141,14 
152,13 
153,14 
152,14 
159,15 
151,15 
160,16 
161,16 
179,17 
171,17 
189,18 
181,18 


199,19 
193 19 


Position 


1828.5181 cm™! 
1828 .5127 
1951.2771 
1906.3849 
1906 .2534 


1843.9477 


1859 .2010 


18 '4.3132 


1889 .3063 


1904.2398 
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